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ABSTRACT 

We report new results obtained from multi-frequency observations of PSR B0826— 34 with 
the Giant Metrewave Radio Telescope (GMRT). (1) We find no evidence of weak emission 
during the typical long null state of this pulsar, simultaneously at 303 and 610 MHz, as well 
as individually at 157, 325, 610 and 1060 MHz at separate epochs. Our limit of non-detection 
is at ~ 1 % or better of the peak of the active state profile, and corresponds to ~ 2 mJy at 
610 MHz. (2) Significant correlation in the total intensity of the individual pulses between 
303 and 610 MHz is reported from the simultaneous dual frequency observations, which is 
indicative of the broadband nature of the emission. We also report correlation between total 
energy in the main pulse and inter-pulse region from the high sensitivity single frequency 
observations at 610 and 1060 MHz. (3) Though we find the drift pattern to be very similar 
in the simultaneous 303 and 610 MHz data, we observe that the drift band separation (P^) 
evolves significantly between these two frequencies, and in a manner opposite to the average 
profile evolution. In addition, we confirm the dependence of Pi on pulse longitude at 303 
MHz and find indications for the same at 610 MHz. We also present results for subpulse 
width (A$ s ) at different frequencies, and as well as a function of pulse longitude. (4) As a 
natural out-come of the simultaneous dual frequency observations, we obtain an accurate DM 
value, equal to 52.2(6) pc/cm 3 , for this pulsar. 

Key words: Stars: neutron - stars: pulsars: general - stars: pulsar: individual: B0826— 34 



1 INTRODUCTION 

The phenomenon of subpulse drifting (first reported by 
Drak e and Craft is manifested as an organised sub- 

pulse behaviour — subpulses appear at progressively changing 
longitude in the pulse window following some particular path. 
The path followed by the subpulses is specific to the individual 
pulsar concerned and is known as drift band. It has been recently 
shown that subpulse drift i ng ma y be fairly common among pulsars 
dWeltevrede etaTI 120061 |2007l) . Wide profile pulsars showing 
this phenomenon can be especially useful tools to understand the 
distribution of the radio emission regions above the neutron star's 
polar cap and the physics of the pulsar emission mechanism. Since 
the wide profiles are usually interpreted as emission from one 
magnetic pole of a highly aligned pulsar (i.e. magnetic axis almost 
parallel to rotation axis), our line of sight samples a large part of 
the polar cap. This leads to enhanced possibility of seeing multiple 
drift bands in such pulsars, which can provide valuable insights 
into the pulsar magnetospheres, as exemplified by some of the new 
and substantive resu lts for PSR B08 26-34 dGupta et al. 1 12004 
lEsamdin et al. 120051). PSR B08 15+09 l lOiao et al. |2004|) and PSR 
B0818-41 tehattacharvva et al. 1120071) . 



ll979l : lBiggs et al. ll985l : lGupta et al. |2004|) have brought out some 
unique properties : strong evolution of the average profile with fre- 
quency, apparent nulling for 70% of time and a remarkable sub- 
pulse drift property — multiple, curved drift bands with frequent 
changes and sign reversals of drift rate. 

Multi-frequency observations, especially when done simulta- 
neously at the different frequencies, are very useful for the study 
of various aspects of the frequency dependence of pulsar radiation 



PSR B0826— 34 is a pulsar with one of the widest known 
pulse profiles. The earlier studies of this pulsar l lDurdin et al. I 



i Taylor. Manchester and Huguenin|[l975l : lBartel and Sieberll978l : 
iBartel Hl98ll : lBhat et al. Ill999l) . Dispersion measure (DM) of pul- 
sars can also be measu red from single ep och simultaneous dual 
frequency observations dAhuia etaL 1 120051) . Simultaneous, multi- 
frequency observations of PSR B0826— 34 can be useful on a vari- 
ety of counts: (i) Determination of accurate DM, as the values re- 
ported and used in the literature span a fairly wide r ange, from 47 to 
65 pel cm 3 . (ii) Study of the pulsar in the null state: lEsamdin et aLl 
d2005l) claim the presence of weak emission, at 1374 MHz, dur- 
ing the long duration nulls of this pulsar. They identify these null 
states as some weak mode of emission where the pulse inten- 
sity is 2% of that of the regular, strong mode, with a pulse pro- 
file that is similar to the strong mode at lower frequencies. Sensi- 
tive, multi-frequency observations at the lower frequencies should 
provide useful information about the behaviour in the null state, 
(iii) Study of simultaneous subpulses: Correlated behaviour of sub- 
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pulse s at different frequencies dTavlor. Manchester and Hugueiifn 



1975 



2003 



Bartel and Sieber 1 119781 : iBartel 1 11981k iKarastergiou et al. 
Kramer etal. I l2003h~ is useful to probe the broadband na- 



ture of the emission mechanism. For pulsars with drifting sub- 
pulses, evolution of P2 (longitude separation be tween adjacent 
drift bands) with frequency can also be signific ant dlzvekova et al. I 
1 199 it iRankin et al. II2005I ; ISmits et all 120071) and a comparison 
of this with the frequency evolution of other profile parameters 
can be useful. Earlier study of drift properties of PSR B08 26— 34 
( iBiggs et al. ll985MGupta et al. ll2004lEsamdin et al. |2005h report 
different P2 values at various individual frequencies which do not 
follow the trend predicted by profile evolution with frequency. In 
addition, for PSR B0826— 34, the intricate nature of the multiple 
drift bands itself deserves a broadband study. 

We observed PSR B0826-34 using the GMRT, simultane- 
ously at 303 and 610 MHz, and individually at 157, 325, 610 and 
1060 MHz in total intensity mode. In this paper we concentrate on 
the results from the simultaneous dual frequency observations, and 
supplement these, where necessary, with the results from the sin- 
gle frequencies observations. In Sect[2]we cover the observations, 
where a new technique for simultaneous dual frequency pulsar ob- 
servations with the GMRT is presented. Different sub-sections of 
Sect[3] describe the data analysis and results for (1) determination 
of DM value from simultaneous dual frequency observations, (2) 
study of the pulsar in its null state and (3) subpulse studies. In 
Sect|4]we discuss the implications of our findings. 



2 OBSERVATIONS 

2.1 Simultaneous dual frequency observations 

The GMRT consists of an array of 30 antennas, each of 45 m diam- 
eter, spread over a region of 25 km diameter, operating at 6 differ - 
ent wave bands between 150 and 1450 MHz dSwarap et al. II 19971) . 
Though designed to function primarily as a aperture synthesis tele- 
scope, the GMRT can be used in an array mode by adding the 
signals from individual dishes, either coherently or incoherently. 
Furthermore, sub-sets of the 30 antennas can be configured com- 
pletely independently in what is referred to as the sub-array mode, 
to effectively provide more than one "single dish", thereby enabling 
pulsars to be observed simultaneously at multiple frequencies. The 
signals from different antennas and observing frequency bands are 
eventually converted to a baseband signal of 16 MHz band width, 
which is then sampled at Nyquist rate. The observing bandwidth 
of 16 MHz is divided into 256 spectral channels in the FX cor- 
relator. These dual polarisation signals from all the antennas are 
then brought to the GMRT Array Combiner (GAC), where they are 
added coherently to get th e phased array outputs for each polari- 
sation dGupta et al. Il200d) . The GAC has the facility of assigning 
different values of gains to each spectral channel of each polarisa- 
tion of each antenna, in a completely independent manner, before 
the addition. 

In the simultaneous multi-frequency sub-array mode of pul- 
sar observations, the baseband signals from antennas in different 
sub-arrays (operating at different radio frequencies) can be added 
together in the same GAC, and the data can be recorded as a sin- 
gle stream of multi-channel data, after integration by the desired 
amount. The dedispersed data streams for the different frequencies 
are then extracted from this single data strea m during off-l i ne pro - 
cessing. Such a technique was employed by lAhuia et al. 1 1 |2005|) . 
who carried out simultaneous dual frequency observations using 



the incoherent array mode of operation of the GMRT. They utilised 
the fact that, while traversing the interstellar medium, different fre- 
quency signals suffer different amount of dispersion delay and as 
a result the pulses arrive at different times (and hence at different 
positions in pulsar rotational phase) at different frequencies. Hence 
the on-pulse data from different frequencies can be recovered from 
appropriate phase windows in the pulse period, during the off-line 
analysis. This scheme provides naturally synchronised data from 
all the sub-arrays, without the need for separate, synchronis ed re- 
ceivers for each sub-array. However, as lAhuia etaL have 
pointed out, this technique works well only for cases where the dis- 
persion curves for different radio frequencies, when mapped to the 
16 MHz baseband, do not overlap with each other for all phases of 
the on-pulse signal. This criterion becomes difficult to meet for pul- 
sars with large duty cycle pulses and for large DMs. Specifically, in 
the case of wide profile pulsars like PSR B0826— 34, this scheme 
will definitely fail due to the fact that the pulses at both frequen- 
cies will overlap in the band and it will be impossible to extract the 
stream of single pulse separately at different frequencies. For such 
pulsars, we have devised a new scheme to carry out simultaneous 
multi-frequency observations. 

As mentioned earlier, the GAC has the facility of assigning 
to each spectral channel, a different value of gain — including 
zero gain. This allows selectable parts of the band to be processed, 
while blocking the rest of the band, for each antenna. The setting 
can be done independently for each antenna. This "band masking" 
technique allows simultaneous multi-frequency observations of all 
types of pulsars, independent of the combination of values of pe- 
riod, pulse width and DM, without having the need for separate, 
synchronised pulsar receiver chain for each sub array. The trick is 
to split the available baseband bandwidth into non-overlapping sub- 
bands, each of these parts being assigned to the antennas of one 
particular sub array, using appropriate band mask settings for each 
antenna. For the antennas of a given sub array, we apply nominal 
gain value for the typical GMRT pulsar observations in the chosen 
part of the band and the rest of the band is masked by applying a 
zero gain value. During simultaneous dual frequency observations, 
the available band is divided into two parts, each part being alloted 
to one of the observing frequencies. 

Using this technique, we observed PSR B0826— 34 simultane- 
ously at 303 and 610 MHz with the GMRT, on October 26, 2003. 
We used two sub arrays, one with 5 antennas (sub array#l), oper- 
ating at 303 MHz, and the other with 10 antennas (sub array#2), 
operating at 610 MHz. We put more antennas at 610 MHz to get 
desired sensitivity, because of the fact that the pulsar is weaker at 
610 than at 303 MHz. The signals coming from the antennas of 
each sub array were added coherently in the GAC. The total band- 
width used for the observations was 16 MHz, which was divided 
equally between the two frequency bands by masking 8—16 MHz 
of the band for the antennas in sub array#l and 0—8 MHz of the 
band for the antennas in sub array#2. Before observing the pulsar, 
both the sub-arrays were individually "phased-up" by estimating 
the phases with respect to a reference antenna in each sub-array 
(from the correlator visibilities recorded on a point source calibra- 
tor) and applying the corrections for these back into the correlator 
hardware. The voltage outputs from the GAC were passed to the 
phased array pulsar receiver where the data were converted to in- 
tensity values and intensities from the two polarisations added to 
obtain the total intensity for each of 256 spectral channels, contain- 
ing information from the two sub-arrays. These data were further 
integrated in time and finally recorded to disk at a sampling rate of 
1.024 ms. During the off-line analysis we further integrated these 
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Figure 1. Upper two panels plot the average profiles (of the active state) for 
PSR B0826-34 at 303 and 610 MHz. The third panel shows the 303 MHz 
profile, aligned to the 610 MHz profile after correction for the dispersion 
delay. MP, IP and OP refer to main pulse, inter-pulse and off pulse, respec- 
tively (as described in the text); MPi and MP2 refer to the two peaks of the 
MP. 



data to a final time resolution of 4.096 ms, as this was adequate to 
study this wide profile pulsar. The duration of the pulsar observa- 
tion was about 1 hr, in which we recorded a data stretch of 2123 
pulses from PSR B0826-34, simultaneously at 303 and 610 MHz. 
In the first 750 of these pulses, the pulsar was in the active state and 
for the remaining time it was in one of its long duration null states. 



2.2 Single frequency observations 

We observed PSR B0826-34 using a large number of GMRT an- 
tennas in phased array mode at 157, 325, 610 and 1060 MHz indi- 
vidually at separate epochs. These observations have significantly 
higher sensitivity (by a factor of 2 or more) compared to the cor- 
responding observations from the dual frequency effort, and were 
carried out primarily to supplement and improve the conclusions 
obtained from the simultaneous dual frequency data set. In addi- 
tion, some of these observations were accompanied with flux cali- 
bration observations using known point sources, which allows for 
absolute calibration of the pulsar's flux. A summary of the main 
parameters of each observing session is given in Table. Q] 



3 DATA ANALYSIS AND RESULTS 

The simultaneous dual frequency data for PSR B0826— 34 were 
dedispersed separately for 303 and 610 MHz, using only the re- 
spective non-masked parts of the baseband signal. The DM value 
used was 52.0, which is close to the final DM value determined by 
us (see Sect l3.lt . The worst kind of radio frequency interference 
that this data is suffering from, is power line related (50 Hz and its 
harmonics). To remove that, the dedispersed data were put through 
a radio frequency interference filtering routine which detected most 
(but probably not all) of the power line interferences and replaced 
them by appropriate random noise. The resulting time series of first 
750 pulses (for which the pulsar was in active state) were then syn- 
chronously folded with the topocentric pulsar period to generate 
the average pulse profiles at the two radio frequency bands, which 
are shown in the top two panels of Fig[T] At both frequencies, the 
emission from the pulsar is present over a wide longitude range, as 
expected. We identify two emission regions of the pulse profiles, 
labeled as main pulse (MP) — roughly from 120°to 300° pulse lon- 
gitude, and inter-pulse (IP) — roughly form 10° to 90°, as marked 
in the panels of Fig[T] In addition, we identify an off pulse region 
(OP) - roughly from 310° to 350° pulse longitude. The MP consists 
of two distinct peaks, denoted by MPi and MP2. The component 
separation between these two peaks is 132° and 1 16° at 303 and 610 
MHz respectively, which is in the same ballpark as reported by ear- 



Una 

lier studies of this pulsar l lBiggs et al. Ill 9851: rGupta et al. 12004]) . 



3.1 Determination of accurate DM 

In order to compute the dispersion delay between the profiles at 
the two frequencies, we need to establish a fiducial point that is 
expected to remain intrinsically fixed at a pulse longitude, even as 
the pulse profile evolves with frequency. For PSR B0826— 34, we 
choose the mid-point of the two peaks of the MP of the average 
profile as the fiducial point. The measured time delay (At) can be 
expressed as, 



At = (N x P) + P E 



(1) 



The topocentric pulsar period P is 1.84881 s. N is the integer num- 
ber of pulsar period delay. Using the observing frequencies, the 
pulsar period and a first guess DM value we estimate TV is equal 
to 1. P t is the delay with in a pulse period. P t is estimated from 
the position difference of the fiducial point between 303 and 610 
MHz pulse profiles. To determine the position of the fiducial points 
at either of the two frequencies, we first estimate the position of the 
two peaks of the MP, by fitting second order polynomial functions. 
Here, the errors in the estimation of the location of the peaks, are 
smaller at 303 MHz, as compared to the 610 MHz, because of more 
complex structure of the profile peaks at 610 MHz (FigO}. Using 
the inferred location of the fiducial points at the two frequencies, 
P e is estimated as 161.1 ms. The error in P t depends on the error 
in the estimation of the fiducial point, which follows from the error 
associated with the estimation of the position of MPi and MP2 at 
the individual frequencies. 

Using Eqn. Q] and the well known formula tha t relates the di ffer- 
ence in travel time at the two radio frequencies dB acker |[l993T) . we 
obtain a value of 52.2 pc/cm 3 for the DM of this pulsar, with an 
error of ± 0.6 pc/cm 3 . The comparison of this result with the DM 
values in the literature is discussed in Sect. 14. II 
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Figure 3. Average profiles for PSR B0826— 34 in the active and null states, from the simultaneous dual-frequency observations, at 303 (left panel) and 610 
MHz (right panel). 
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Figure 4. Same as Fig.[5]for higher sensitivity single frequency observations at 325 (left panel) and 610 MHz (right panel) 
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Figure 5. Same as Fig.|3]for higher sensitivity single frequency observations at 157 (left panel) and 1060 MHz (right panel) 
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Table 1. Table containing, (a) summary of all the observations presented in this paper, (b) rms values for pulsar active and null state, (c) the flux of the peak of 
the pulsar active state and the non-detection limit for null state and (d) the mean flux for the FP, MP and IP regions of pulsar active state for all the observing 
sessions (see Fig.[3]for the positions of the MP, OP and IP windows). 



Frequency 


State 


No of pulses 


No of 


rms 


• 


Flux of peak 


Non-detection 




Mean Flux 




(MHz) 


(Active/Null) 


analysed 


antennas 








limit 














used 


(mJy) 


(mJy) 


(mJy) 




(mJy) 












MP 


OP 






FP 


MP® 


IP® 


303+ 


Active 


750 


5 


t7 1 


28 


1687* 


- 


443* 


754 * 


99* 




Null 


750 


5 


28 


28 


" 


3 1 * 








610t 


Active 


750 


10 


S7 


j 


250* 


- 


119* 


178* 


123* 




Null 


750 


10 


4 


4 




5* 








157* 


Active 


970 


16 


427 


1 3 


1780* 


- 


370* 


656* 


42* 




Null 


970 


16 


13 


13 




23" 








325 § 


Active 


940 


21 


304 


5 


1330* 


- 


393* 


763* 


62* 




Null 


487 


17 


5 


5 












610* 


Active 


705 


17 


85 


2 


410*° 




156*° 


248*° 


166*° 




Null 


705 


17 


2 


2 












1060* 


Active 


2900 


21 


120 


10 


1310* 




479* 


485* 


989* 




Null 


2900 


21 


10 


10 




15* 









f : Active and null state data recorded from same simultaneous dual frequency observing session 
| : Active and null state data recorded from same single frequency observing session 
§ : Active and null data recorded from single frequency observations from different days 

• : Value quoted for profile data with time resolution of 4.096 ms 

o : 3-sigma Non-detection limit for the null state, for profile data with time resolution of 20.048 ms 

* : Flux estimated using standard GMRT system parameters using relative calibration procedure 
o : Flux estimated using absolute calibration with a flux calibrator source 

: Mean flux is calculated up to 50% intensity lavel for MP and IP 



3.2 Study of the pulsar in null state 

Both short (~ few pulses) and long duration (~ few hundred pulses 
or more) nulls are present in our simultaneous dual frequency data. 
FigH] shows the subpulse patterns at 303 and 610 MHz. For the 
initial 750 pulses, the pulsar was mostly in an active state except 
for a few short duration nulls which appear to be simultaneous at 
the two frequencies (e.g. pulse # 250 to pulse # 260, pulse # 414 to 
pulse #421, pulse #613 to pulse #621 and pulse # 683 to pulse # 
691, in FigHJ. 

After the first 750 pulses, the pulsar goes into a long, appar- 
ently null state, and this happens simultaneously at both the fre- 
quencies. The pulsar remains in this state for the rest of the obser- 
vations (about 1373 pulses). We obtain the average profile of the 
active state by integrating the first 750 pulses. Similarly, an average 
profile of the apparent null state of comparable duration is obtained 
from the next 750 pulses (see Fig|3j. We see that during the null 
state, there is no evidence of any emission that matches with the 
active state average profile for the MP or IP regions, at either of the 



two frequencies. The mean power level in the null state profiles is 
comparable to the mean level of the OP region of the corresponding 
active state profiles. The fluctuations of the signal level in the null 
state profiles, though not completely Gaussian noise-like (proba- 
bly due to low-level unexcised interference and other sources of 
"red" noise), show no correspondence with the features in the active 
state profiles. Furthermore, these features in the null state profiles 
have no resemblance or correlation between the simultaneous data 
at 303 and 610 MHz, indicating that these are not genuine profile 
features. In order to further check for any short duration bursts of 
strong emission in the null state, we divide the data (the first ~ 1500 
pulses in the observation) into 20 blocks and compared the average 
profiles from these blocks. For the first 10 blocks, during which the 
pulsar is in the active state, we see the characteristic pulse profile at 
either of the frequencies. For the next 10 blocks, during which the 
pulsar is in the null state, we do not see any systematic signature 
of emission, at either of the two frequencies, in any one of the 10 
blocks. 

For a more thorough investigation of the null state, we com- 
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pare the rms values for the MP, OP and full pulse (FP) windows 
of the integrated pulse profiles (see Table. [TJ. The rms during the 
pulsar null state is similar for all the windows (MP, OP and FP) and 
is also similar to the rms of the OP region of the pulsar active state. 
These rms values are much smaller than the rms of the MP window 
of the active state, and their ratios to the corresponding mean val- 
ues are, to first order, consistent with thermal noise statistics. This 
behaviour is consistent for both the frequencies, indicating that the 
null state signal is similar to the OP region of the pulsar's active 
state, and contains no detectable level of emission from the pulsar. 
To further constrain the above conclusion, we integrated the pro- 
file data to a time resolution of 20.048 ms, which is a reasonable 
compromise between integrating all of the pulsar's signal into a 
few bins and retaining enough time resolution to detect any large 
scale emission structure. Comparing the peak deflection in the ac- 
tive state profiles with three times rms of the fluctuations in the null 
state profiles, we conclude that there is no emission in the null state 
down to a level of ~ 2% of the peak in the active state, for either 
frequency. 

The dual frequency observations used a limited number of an- 
tennas at each frequency and are not very sensitive. In order to 
achieve better sensitivity, we use the results from the single fre- 
quency observations carried out with larger number of antennas 
in phased array mode. Fig. [4] and Fig. [5] show the active and null 
state profiles from these observations. Whereas the active and null 
state profiles for the 325 MHz observations are from different days, 
those for the 157, 610 and 1060 MHz observations are cases where 
the pulsar transited from the active to null state during the same 
observing session. 

As can be seen, there is no evidence for systematic emission 
from the pulsar in the null state at any of the frequencies. Non- 
random features in the null state profiles, where seen, have no cor- 
respondence between different frequencies, or even for the same 
frequency at different epochs. 

Comparing the peak deflection in the active state profile with 
three times the rms of the fluctuations in the null state average pro- 
file at a time resolution of 20.048 ms, we conclude that at 610 MHz 
there is no emission in the null state down to a level of ~ 0.5 % of 
the peak in the active state. At 1060 MHz no signature of emission 
is seen in the null state profile down to ~ 1.2% of the peak in the 
active state. Similarly, at 157 MHz the limit is ~ 1.3%. 

From the absolute flux calibration available, we determine the 
flux of the peak for 610 MHz active state to be ~ 410 mJy, and the 
3-sigma non-detection limit in the null state corresponds to about 
2 mJy. Similarly, from the absolute flux calibration available for 
the null state observations at 325 MHz, we find the 3-sigma non- 
detection limit at this frequency to be ~ 7 mJy. For the cases where 
we do not have any flux calibrator source observed, we follow a rel- 
ative calibration procedure, using the knowledge of the observing 
parameters and the background sky temperature near the source. 
As a crosscheck we find that, for the single frequency 610 MHz 
observations, the flux values estimated with absolute and relative 
calibration are similar. 

Though the scintillation properties of PSR B0826— 34 is not 
studied, still we try to understand the effect of scintillation on the 
observed flux, by comparing with two nearby pulsars for which the 
scintillation properties are known. The mean flux densities for two 
nearby pulsars, PSR B0823+26 and PSR B0834+06, a re observed 
to ch ange by ~ 35% between the observing epochs dBhat et al~l 
Between the two observing sessions of PSR B0826— 34 at 
610 MHz, we observe, 27% change of the mean flux of the FP. 
The mean flux at 610 MHz is close to the mean flux estimated by 
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Figure 6. Spectra from calculated mean flux densities (Table. [T} for FP 
(solid), MP (dashed), IP (dotted) 

lBiggsetal.Ul985h at 645 MHz. Estimated mean flux changes by 
12% between two observing sessions at 303 and 325 MHz. So, the 
flux variation observed by us are well within the 3 a limit of what is 
reported for the nearby pulsars. Hence, it is likely that the estimated 
mean flux values of PSR B0826— 34 are affected by the scintillation 
to a moderate extent. 

In Table[T] we list the flux at the peak and the mean flux of 
PSR B0826-34 for FP, MP and IP regions at different frequencies. 
Fig.|6]illustrates the frequency evolution of the mean flux density of 
the FP, MP and IP regions. At lower frequencies like 157 MHz, the 
IP is almost absent and the MP is the dominant component. With 
increasing frequency the mean flux density of the IP increases. At 
610 MHz, the mean flux density of the MP and the IP are almost 
equal and at 1060 MHz, the IP flux density dominates the MP. For 
the simultaneous dual frequency observations, we see that the FP 
and MP has quite steep spectral index, a/ p ~ -1.9 and a m p~ -2.1 
,whereas the IP has flatter spectral index, c*i P ~ 0.3 (where a is 
defined by I oc u a ). These are similar to the results reported by 
iBiggs et al~1 jl985h between 408 and 645 MHz. From single fre- 
quency observations at 610 and 1060 MHz, the calculated spectral 
indexes are, 2, a mp ~ 1.2 and ai p ~ 3.2. Hence, the spectral 

index a for this pulsar is different for different profile components 
and also evolves with frequency. This can be understood by con- 
sidering that the IP emission is coming from a second ring of emis- 
sion which is missed by our line of sight at the lowe r frequencies 
and is visible at the higher frequencies, as proposed bv lGupta et al. I 
d2004h . 

To summarize, our results show similar limits (~ 1% of peak 
of active state profile, or better) of non-detection of emission in 
the null state, at frequencies of 157, 303, 325, 610 and 1060 MHz. 
These results are in apparent disagr eement with those obtained at 
1374 MHz bv lEsamdin et al. I d2005l) , who report presence of weak 
emission in the apparent null state. We discuss the implications of 
this in Sect|4~2l 

3.3 Study of subpulse emission 

We align the 303 and 610 MHz data from the simultaneous dual 
frequency observations, by removing the initial 2009 ms of data 
(corresponding to the dispersion delay) from the 303 MHz data. 
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Pulse Phase (deg) Pulse Phase (deg) 



Figure 8. Gray scale plot of single pulse data for pulse number 400 to 600, from the simultaneous dual-frequency observations, at 303 (left panel), and 6 1 
MHz (right panel). 




Lag (Pulsar Period) 

Figure 7. Plot of correlation function versus pulse lag: autocorrelation func- 
tion of TE mp for 303 MHz (solid), same for 610 MHz (dashed), autocorre- 
lation function of TEjj, at 610 MHz (dotted), crosscorrelation functions of 
the TE mp between 303 and 610 MHz (dash-dot). The plot for the autocor- 
relation function starts from unit lag for all the three cases. 



This allows us to directly compare the corresponding single pulses 
from these two frequencies. 



3.3.1 Pulse energy correlation 

To study the correlation of pulse energy fluctuations as a function 
of frequency and pulse lag, we compute the correlation of the total 



energy under the MP (hereafter TE mp ), for different pulse lags, 
individually at 303 and 610 MHz, as well as between these two 
frequencies. For a train of single pulses with pulse numbe r k, 1 < 
k < M, the correlation function is given by dBartelll98lh . 



C(0 



(a 2 h a 



1/2 



(2) 



where i, j denote the concerned frequencies (303 or 610 MHz); 
If i (k) and //. (k) are the TE mp of the k th pulse for frequencies 
i and j, respectively; (// 4 ), (//■) are the mean energies, and Of t , 
o~f. are the variances for TE mp over M pulses. A similar procedure 
is adopted for the IP, where applicable. 

Fig[7J shows the correlation function (for pulse number 1 to 
750), for autocorrelation of TE mp at 303 and 610 MHz, autocorre- 
lation of TEip at 610 MHz and crosscorrelation of TE mp between 
303 and 610 MHz, as a function of increasing pulse lag. The TE mp 
for the active state (pulse number 1 to 750), with zero pulse lag 
(i.e. I = 0), is 47% correlated between 303 and 610 MHz. More- 
over, this level is very similar to the correlation at each of the two 
frequencies. All these three curves (as well as that for the autocorre- 
lation of the IP at 610 MHz), follow very similar shapes with pulse 
lag, and show positive correlations out to fairly large lags: about 
1 1 periods. Similar results are obtained from correlation studies of 
the higher sensitivity single frequency data. Furthermore, from the 
single frequency data, we find a crosscorrelation between TE mp 
and TE ip of about 22% at 610 MHz and 17% at 1060 MHz. In 
the autocorrelation function for 303 MHz (Fig. [7), and also in the 
autocorrelation of the 610 MHz data from the single frequency ob- 
servations (not shown in Fig. [7} a secondary maximum is present 
at pulse lag / = 13, which is not consistently seen in all correlation 
functions. These different aspects of the correlation results, includ- 
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ing comparison with existing results for other pulsars, are covered 
in SectgJ 



3.3.2 Study of subpulse drifting 

Drifting subpulses observed simultaneously at 303 and 610 MHz 
(e.g. Fig|2j, broadly reproduce the features reported by ear- 
lier studies at various individual frequ encies dBiggs et al. 1 1 19851 : 
IGupta et al. 2004; Es amdin et al. 120051) . For a closer investigation 
of the simultaneous subpulse behaviour, we zoom into a sequence 
of 200 pulses (pulse # 400 to pulse # 600) as shown in Fig[8] At 
any given time, the simultaneous multiple subpulses present in the 
main pulse window follow the same drift rate and sign, at both the 
frequencies. The pulse regions showing different drift rates of op- 
posite signs are always connected by a region showing smooth tran- 
sition of drift rate, as reported bv lGupta et al. I ( 12004 ). Such transi- 
tions, seen in pulse numbers 440 to 460, 500 to 520, 540 to 560 
and 570 to 590, also occur simultaneously at 303 and 610 MHz. 
The gray scale plot of the single pulses from the higher sensitiv- 
ity 610 MHz single frequency observations (Fig.[9]l show coherent 
drifting in the MP and IP regions, with approximately 6 drift bands 
in the MP and 4 drift bands in the IP. Subpulse drifting is observed 
under the two peaks of the MP at 157 MHz (Fig. [9}. Though the 
drifting subpulses show very similar structure in our simultaneous 
dual frequency observations, there are subtle differences such as 
changes in the separations and widths of the subpulse as a function 
of frequency. To investigate the change of drift band separation, we 
compute P2 from the autocorrelation function of the single pulses, 
averaged over the 750 pulses in the active state. In order to clearly 
see the secondary maximum in the correlation function (which will 
give us the P2 value), we correlate adjacent pulses with different 
longitude lags, rather than the same pulse (which is the traditionally 




Lag (deg) 

Figure 10. Autocorrelation results for section of profile for the MP (longi- 
tude 120° to 300°), calculated with one pulse offset, for pulse # 1 to 750 at 
303 (solid) and at 610 MHz (dash). Secondary maximum of the autocorre- 
lation function for 303 MHz is 22.6°and 610 MHz is 27.9°. 



correct procedure). This one pulse offset method suppresses a large 
amount of the correlation within the same subpulse, and as a conse- 
quence the secondary peak becomes stro nger and easily detec table. 
This can be seen by comparing Fig.4 of IGupta et al. I d2004l) with 
our Fig[l0] We find P 2 = 22.6±0.8°at 303 MHz and 27.9±0.8°at 
610 MHz. Thus, the value of P2 evolves with frequency. Any doubt 
about epoch dependent variations of P2 manifesting as variation of 
P2 with frequency can clearly be ruled out with these simultane- 
ous dual frequency results. Further frequency evolution of P2 is 
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Table 2. DM values of PSR B0826-34 from the literature survey. 



Reference 



DM value quoted 

(pc/cm 3 ) 



Durdinetal. (19791 
EPN database 

Gould and Lvne (19981 & D'Amico et al. (1998) 



Taylor, Manchester and Lyne ( 1 993) 
Gupta etal. (2004) 
Hobbsetal. (2004) 
Esamdinetal. (2005) 
This paper 



S2t 
47+ 

47±4* 
47t 
52.9±0.6* 

65.6t 
52.2±0.6§ 



f : DM determination method not described 
\ : DM determined using timing analysis 

§ : DM determined using Simultaneous dual frequency observations 



quency, we compute the autocorrelation function (without any 
pulse offset) and find the full width at half maximum (FWHM). Av- 
erage value for A$ s is found to be 17.4° at 303 MHz, and 14.2° at 
610 MHz, again confirming a frequency evolution. We collect val- 
ues of A$ s at different frequencies, from our observations as well 
as from literature, in Table. [5] We observe longitude evolution of 
A$ s : for region I (from 118°to 197°pulse longitude) A$ s is sig- 
nificantly more than that for region II (from 197° to 276° pulse lon- 
gitude) for all the frequencies. For a more detailed study of the 
longitude evolution of A& 3 within the MP, we followed a window 
scheme, similar to that employed for the P2 variation study. We 
find, A$ s varies smoothly over the pulse window, in a manner that 
is somewhat similar to that of Pi variation with longitude, lower 
A$ s values in the middle part of the MP and higher values to- 
wards the edges (A<J> 3 values near MPi being more than those near 
MP 2 ). 



obtained from our single frequency observations at different fre- 
quencies, the results from which are summarised in Table. [4] In 
addition, we find P2 ~ 20° for the IP region, from 610 and 1060 
MHz single frequency observations. 

To study the longitude evolution of P2, we repeat the correla- 
tion analysis over few chosen windows (~ 32° wide) of pulse lon- 
gitude which are shifted systematically over a selected portion of 
main pulse window, in steps of 8° . Fig[TT]is a sample plot of corre- 
lation functions for 7 different pulse windows at 303 MHz and 610 
MHz. For 303 MHz we observe a smooth trend for variation of P2 
across the pulse window, with a minimum of 20.4° to a maximum of 
about 25°, which is compatible with the mean of 22.6°we obtained 
from the full main pulse. P2 increases from the centr e towards the 
edge o f the MP, similar to the variation reported by iGupta et al. I 
d2004l) at 318 MHz. For 610 MHz, the variation is from a mini- 
mum of 26.3°to a maximum of 30.3°, again with a tendency to be 
lower in the middle of the MP, though the results are not as clear as 
for 303 MHz. 

To study the variation of subpulse width, A$ s , with fre- 



4 DISCUSSIONS AND SUMMARY 

4.1 Determination of accurate DM value 

Table. [2] compares our result for the DM value of PSR B0826-34 
with other available values from the literatur e — these cover a 
wide range, with the two most recent papers <Guptaetal.ll2004 
lEsamdin et al. I l2005t) using values of 47 and 65.6 pc/cm 3 . It is 
clear that, unlike most normal pulsars, DM determination for this 
pulsar is a difficult and trick exercise, mainly because the profile 
is quite complex, very wide and strongly evolving with frequency. 
Given our more accurate method of determining the DM, we be- 
lieve that the either of the extreme values above are unlikely to be 
co rrect. We note that our DM result is closest to value obtained 
bv lHobbs et all | |2004|) . using multi epoch timing. The accuracy of 
both these methods depends on the choice of fiducial point. The 
advantage of our method of DM determination is that, the obser- 
vations at a single epoch are self sufficient for obtaining the DM 
value at that epoch, without requiring any absolute timing infor- 
mation. This method can facilitate the DM determination of other 
complex and wide profile pulsars. 
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Table 3. Correlation of pulse energies between different frequencies of observations. 



Pulsar 


Frequency range 
MHz 


Correlation Coefficient 


Pulse Lag up to which 
correlation is positive 




Reference 




B0031-07 


275 


70% 




Taylor, 


Manchester and Hueuenin 


(1975) 




415 












B0031-07 


102 
406 


50% 


iot 




Izvekova et al. (1991) 




B0329+54 


327 
2695 


75% 


2t 




Battel and Sieber 0978) 




B0329+54 


102 
1720 


50+30% 


4t 




Bartel (1981) 




B0329+54* 


238 
626 


67±3% 






Kramer et al. (2003) 




B0809+74 


275 


70% 




Taylor. 


Manchester and Hueuenin 


(1975) 




415 












B0809+74 


102 
1720 


30±10% 


4t 




Bartel (1981) 




B0826-34 


303 
610 


47% 


11 




This paper 




B0834+06 


275 


70% 




Taylor, 


Manchester and Hueuenin 


(1975) 




415 












B1133+16 


275 


78% 




Taylor, 


Manchester and Hueuenin 


(1975) 




415 












B1133+16 


327 
2695 


75% 


2t 




Battel and Sieber (1978) 




B1133+16§ 


341 
626 


72+4% 






Kramer et al. (2003) 




B 1508+55 


102 
1720 


30±20% 


It 




Bartel (1981) 





- : indicate that value for the corresponding parameter is not known 

f : values obt ained from visual ins pection of plots for crosscorrelation coefficients vs pulse lag from respective references 
% : Table 4 oflKram er et al.l fe003l) list the crosscorrelation coefficients between various other frequencies for B0329+54. 
Kara stereiou et al. I i2003l) reported similar values of crosscorrelation coefficients 

§ : Table 5 o f lKrameretaU 120031) list the crosscorrelation coefficients between various other frequencies for Bl 133+16 



4.2 Study of the pulsar in null state 



We see short duration nulls of a few pulses, which appear to be 
simultaneous over the frequency range of 303 to 610 MHz. In ad- 
dition, this pulsar exhibits long durations of what is apparently a 
null state, which also appears to be broad band over this frequency 
range. However, we find no evidence for any weak emission in 
these long, apparently null states of the pulsar. Non-detection of 
any correlated emission structure in the null state profiles from the 
simultaneous dual frequency observations (at 303 and 610 MHz) 
are a clear proof of this. For frequencies above (up to 1060 MHz) 
and below (down to 157 MHz) also, we see no evidence for any 
weak emission. Our best limits on this non-detection are compa- 
rabl e to, or a few times b etter, than the level of detection (2%) 
that lEsamdin et aL 1 d2005h report for the null state of this pulsar, 
at 1374 MHz. These comparisons are all with respect to the peak of 
the emission profile in the active state. We have also obtained abso- 
lute flux limits for the non-detection at various frequencies, which 



should be a useful comparison standard for any more sensitive stud- 
ies in the future. 

What could b e the reason for the di fference between our re- 
sults and those of lEsamdin et al. ] d2005h ? It could be a case of 
frequency dependent nulling, as reported bv lBhat et al. I d 19991) for 
PSR B 1133+ 16 where the pulsar exhibits cases of short duration 
nulls at lower frequencies like 325, 610 and 1400 MHz, while there 
is visible emission at 4850 MHz. However we believe that this ex- 
planation is difficult, given that we see no evidence for emission till 
the frequency of 1060 MHz, coupled with the fact that the emission 
properties of this pulsar appear to be significantly correlated over 
an octave range in freque ncy. Hence, we are una ble to reconcile our 
results with the results of lEsamdin et al. I d2005l) . 

4.3 Study of subpulse emission 

(i) Pulse energy correlations : Our result that the fluctuations of 
the main pulse energies at 303 and 610 MHz are correlated up to 
47% level, indicates that the emission process generating the single 
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Table 4. Frequency dependence of profile width and P 2 for different pulsars. 



Pulsar 


Frequency 


Profile width 


P 2 


Comments 


Reference 






A$° 










MHz 


(°) 


(°) 






B0031-07 


60 


75 


20.6 


P 2 oc y-°- 05 


Izvekova et al. (1991) 




102 


65 


19.8 


A$ oc u~ 0A 




406 


40 


18.7 







B0031-07 157 44+5 18.9±1.2 P 2 evolution Smits et al. (2007) 

A-mode 



157 


44+5 


18.9±1.2 


P 2 evolution 


243 


39.5+1.5 


24±4 


does not follow 


325 


39.4±0.3 


19.8±1.0 


definite trend 


607 


32±4 


21+3 




840 


32.2±0.8 


10+3 




1167 


25±6 


14+5 




4850 


23.6+0.6 


14+2 





B0031-07 


157 


46.3+1.4 


19.2+0.7 


P 2 decreases 


B-mode 


243 


42.3+0.6 


18.7+1.6 


or remains unchanged 




325 


40.3+0.12 


18.7+0.6 


with increasing 




607 


39.3+0.9 


18.7+0.6 


frequency 




840 


43+0.9 


14+1.6 






1167 


38+5 


12+3 





Smits etal. (2007) 



B0320+39 


102 


11 


3.8 


P 2 oc v-O-iB 


Izvekova et al. (1991) 




406 


8.5 


3.1 


A$ oc u-° A 


B0809+74 


102 


34.8 


14.7 


see Sectl4~3l 


Daviesetal. (1984) 




406 


26.0 


10.8 








1412 


29.8 


8.6 








1720 


21.2 b 


8.1 




Bartel (1981) 


B0826-34 


157 


140+4+ 


23 ±0.8 




This Paper 




303 


132±4t 


22.6+0.8 


A* oc v~ - 2 


This paper 




318 


134±4t 


24.9+0.8 


P2 increases with 


Guntaetal. (2004) 




610 


116+4+ 


27.9+0.8 


increasing frequency from 


This paper 




645 


112+4+ 


29+2.0 


303 to 645 MHz 


Biggs etal. (1985) 




1060 


117+4+ 


28.7+0.8 




This Paper 




1374 




27.5 




Esamdinetal. (2005) 


B 1237+25 


430 


16 a 


10.6+2 


P 2 oc v -° 16 


Wolszczan et al. (1981) 




1700 


14.5 6 


9.1 + 1.5 


A* oc iz-0.07 





B2016+28 


102 


20 


S.5 


P 2 oc u~ 01 


Izvekova et al. (1991) 




406 


15 


6.2 


A$ oc v~ - 2 






1412 


15 


5.5 






B2020+28 


430 




28.3+4 


P2 increases with 


Wolszczan et al. (1981) 




1700 


I6.2 b 


51.4+7 


frequency 





o : Profile width at 10% intensity level 

a : average profile width at 400 MHz from iGould and Lvnel 1998) 



b : average profile width at 1600 MHz from iGould and Lvne Hl998 | ) 
1,11 : left and right part of MP respectively, see Fig. 1 of lBiggsetal 



1985) 



f : separation between two components of main pulse 



pulses has a high degree of correlation over a this octave bandwidth. 
What is equally noteworthy is that this level of correlation is pretty 
much the same as that at each of these individual frequencies, indi- 
cating that there is practically no loss of correlation over this range 
of frequency. Table[3] lists the correlation coefficient of pulse ener- 
gies for different pulsars between two different observing frequen- 
cies, reported in the literature. All the pulsars listed there show sig- 
nificantly high correlation, supporting the case for the broad band 
nature of pulsar radio emission. At the same time, the significantly 
lower correlation between the main pulse and inter-pulse energies 
at 610 and 1060 MHz, implies that the physical processes respon- 



sible for the MP and IP emission may be much more unrelated to 
each other. 

The result that the intensity correlations are positive for large 
lags implies that there is some kind of memory in the underly- 
ing subpulse structure. That this memory is the longest for PSR 
B0826— 34 (amongst all known cases), is quantified in column 4 of 
Table(3] This is likely to be related to the large number of subpulses 
(or drift bands) that each main pulse energy estimate encompasses. 
The memory could then correspond to the time taken for the emis- 
sion from a given subpulse spark column to drift through the main 
pulse window. 



© 2007 RAS, MNRAS 000.QJQ4] 



Results from multi-frequency observations ofPSR B0826—34 13 



The other memory that the main pulse intensity variations ex- 
hibit is the intriguing feature of the peak at lag 13 seen in some 
of the correlation results. This could be an indication of the total 
time for rotation, around the magnetic axis, of the pattern of drift- 
ing subpulses — the caro usel rotation time. F or example, in the 14 
spark model presented bv lGupta et al. the drift rate is such 

that one subpulse drifts to a location close to that of the adjacent 
subpulse in one pulsar rotation, leading to a carousel rotation time 
of Pa = 14Pi. However, if we believe that there are 13 sparks, 
as indicated in the model given bv lEsamdin et al. |i 20051), w ith the 



Table 5. Frequency dependence of subpulse width (A$ s ) of PSR 
B0826-34. 



1 20041) , then 



same drift rate condition as proposed bv lGupta et al. 
we would have P4 = 13Pi, matching with the secondary peak 
feature seen in our data. 

(ii) Study of subpulse drifting: PSR B0826-34 exhibits sub- 
stantial variations of P2 at different frequencies, usually reported at 
different observing epochs (Table The different values of P2 de- 
termined from the simultaneous data at 303 and 610 MHz strongly 
supports the claim that these are genuine variations with frequency, 
and not due to some time dependent phenomenon. The collected 
results in Table|4]are consistent with a scenario where P2 increases 
with frequency from 157 to 645 MHz, but then appears to satu- 
rate or even turn over at the higher frequencies like 1060 and 1374 
MHz. For the low frequency range, we find that this frequency evo- 
lution (P2 oc v 3 ) is of opposite sense to that for the profile width, 
as measured by the separation between two components of the MP 
(A<£> oc v~ ' 2 ). Only two other pulsars in Table|4]show such an ab- 
normal behaviour : PSR B2020+28 shows a clear opposite sense 
(though only 2 frequency points are available) and PSR B0031— 07 
in the A-mode shows erratic increases and decreases of P2 with 
increasing frequency. 

For most pulsars, though both P2 and the average profile 
width decrease with increasing frequency, the evolution of these 
two q uantities with frequenc y does not exactly match (see Ta- 
ble|4ll. Ilzvekova et all J 1 99 lh showed in detail that P2 is, in gen- 
eral, less dependent on frequency, as compared to the profile width 
— this is reflected in our Table|4]as well. T his has been explained 
l lGil and Krawczvk II 19961 : iGil et al. Il2002h to be a natural conse- 
quence of the conal model of pulsar beams l lRankin et al. II 1993h : 
the subpulse enhancement follows a narrow bundle of dipole field 
lines, while the mean enhancement corresponding to the peaks of 
mean profile components, is distributed over a cone (or a number 
of nested cones) of dipole field lines. Thus, the frequency depen- 
dence of subpulses and pulse widths should in general be different. 
However, these models do not appear to explain the opposite sense 
that we see in PSR B0826-34. 

Another pulsar with anomalous frequency variations of pro- 
file width and P 2 is PSR B0809+74. P 2 being 1:8 times greater 
at 102.5 MHz t han a t 1720 MHz, for this pulsar, is explained by 
lEdwards et al. U2003I) as an artifact of t he phase step, which is only 
present at high frequencies. Whereas I Rankin et al. ] l l2005l) anal- 
ysed the P2 variations between 40 to 1400 MHz for this pulsar, 
and showed that they arise from incoherent superposition of two 
orthogonal modes of polarisation. It will be interesting to see if de- 
tailed polarisation studies of PSR B0826— 34 can throw some light 
on the P2 variations with frequency. 

The other interesting result is the variation of P2 with 
pul se longitude, where we confirm the behaviour first reported 
by iGupta etal.l d2004l) at 318 MHz (a variation from 21.5°to 
27° within the MP), and a lso find some signs f or a similar ef- 
fect at 610 MHz. Similarly. lEsamdin et all d2005h re ported that at 
1374 MHz, P 2 varies from 26.8°to 28°within the MP. lGupta etaTl 



Frequency 
(MHz) 


MP 


A<3? s 

(°) 

I 


11 


IP 


Reference 


157* 


15.8 § 


22. 1 § 


10.2 § 




This paper 


303+ 


17.4§ 


28.4§ 


12.6 s 




This paper 


610+ 


14.2§ 


22§ 


12.6 § 


12.6 § 


This paper 


610* 


14.2§ 


22.9 § 


12.6 § 


12.6 § 


This paper 


645 




14* 


5.9* 


4 * 


Bisesetal. (1985) 


1060* 


14.2§ 


20.5 § 


12.6 § 


17.4§ 


This paper 


1374 


13.1±1.9 






19±3.5 


Esamdinetal. C2005) 



f : From simultaneous dual frequency observing session 
\ : From single frequency observing session 
1 : 1 18°to 197°pulse longitude 
II : 197°to 276°pulse longitude 

* : Lag corresponding to half width at half maximum (HWHM) values 
5 : Error associated in determination of A<l> s is ±1.6° 



d2004l) have provided an explanation of this in terms of the geomet- 
rical effect of the way our line-of-sight traverses the polar cap. 

A<E>s determined in this work is compa red to the result quoted 
by dBiggs et al. Hl985l : lEsamdin et al. 1120051) at 645 and 1374 MHz 
respectively in Table. [5] AQ S is observed to be frequency depen- 
dent, average A$ a at 303 MHz is more than 610 MHz value (see 
Sect l3.3t . In addition, A$ s is observed to be longitude dependent: 
A«I> S in region I is more than region II (Table. f5). Similar longitude 
evo lution of A$ a within t he MP is reported bv lBiggs et al. 1 1 19851) 
and lEsamdin et al. I d2005l) . 

To summarise, using results from high sensitivity single pulse 
observations simultaneously at 303 and 610 MHz, and individu- 
ally at 157, 325, 610 and 1060 MHz, we have investigated different 
aspects of PSR B0826— 34: (1) precise value of DM (2) emission 
properties of the long duration null states and (3) simultaneous sub- 
pulses and details of the drifting behaviour. 
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